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ABSTRACT 

The  effects  of  variable  properties  on  heat  transfer  in  a 
film  cooled  turbulent  boundary  layer  were  investigated.  A 
new  procedure  was  developed  to  deduce  adiabatic  effective- 
ness from  heat  transfer  coefficients  based  on  the  wall  to 
freestream  temperature  difference,  where  both  are  represen- 
tative of  variable  property  flow  conditions.  The  new  tech- 
nique was  shown  to  be  valid  using  data  from  the  literature 
for  injection  into  a  turbulent  boundary  layer  from  one  and 
two  rows  of  injection  holes.  From  these  results,  the  varia- 
tion of  the  coolant  to  mainstream  density  ratio  was  shown  to 
have  a  significantly  greater  effect  on  heat  transfer  than 
variations  of  viscosity  and  thermal  conductivity. 
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NOMENCLATURE 

d      -  Injection  Hole  Diameter 

H     -  Enthalpy 

h     -  Heat  Transfer  Coefficient  with  Film  Cooling 

(spanwise  averaged),  q"/ (T^^  -  T^)  and  q"/(T^^  -  T^^J 

h£    -  Heat  Transfer  Coefficient  with  Film  Cooling 
(spanwise  averaged),  q"/ (T^^^.  -  T^) 

hg     -  Heat  Transfer  Coefficient  without  Film  Cooling 

(spanwise  averaged),  q'V(T^^  -  T^)  and  q"/ (T^^  -  T^) 

2       2 

I  -  Injection  Momentum  Flux  Ratio,  P^U^/p^U^ 

k  -  Thermal  Conductivity 

m  -  Mole  Fractions  of  Constitutant  Gases 

M  -  Blowing  Ratio,  P^U^/p^U^ 

P  -  Static  Pressure 

Pr  -  Molecular  Prandtl  Number 

q"  -  Wall  Heat  Flux 

Re  -  Reynolds  Number,  Pc^'^c^^'^c 

T  -  Absolute  Temperature 

U  -  Mean  Velocity 

u  -  Velocity  Component  in  X  Direction 

V  -  Velocity  Component  in  Y  Direction 

W  -  Molecular  Weights  of  Constitutant  Gases 

X  -  Distance  from  Downstream  Edge  of  Injection  Holes 
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x/d    -  Nondimensional  Distance 

a     -  Thermal  Diffusivity 

p     -  Density 

n     -  Adiabatic  Film  Cooling  Effectiveness  ( spanwise 
averaged) ,  eqn.  15  and  22 

y     -  Absolute  Viscosity 

V  -  Kinematic  Viscosity 

0     -  Nondimensional  Coolant  Temperature,  eqn.  16 

Y  -  Injection  Angle  Measured  from  Mainstream 

Flow  Direction 

C      -  Nondimensional  Distance  Parameter,  eqn.  25 
X  -  Second  Coefficient  of  Viscosity 

SUBSCRIPTS 

aw  -  Adiabatic  Wall 

c  -  Coolant 

e  -  Some  Reference  Value 

o  -  Total  Condition 

r  -  Recovery  Condition 

°*  -  Freestream 

SUPERSCRIPTS 
—      -  Averaged  or  Mean  Value 

-  Variable  Property 
*      -  Fluctuation  from  the  Averaged  Value 
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I.  INTRODUCTION 

Over  the  past  thirty  years,  gas  turbines  have  become 
smaller  and  lighter  with  improved  fuel  economy.  These 
desirable  characteristics  have  been  achieved  by  raising  the 
turbine  entry  temperatures.  To  prolong  useful  blade  life  in 
this  high  temperature  enviornment,  cooling  systems,  thermal 
barriers,  and/or  different  materials  capable  of  withstanding 
higher  temperatures  are  utilizied. 

This  thesis  concentrates  on  film  cooling.  Film  cooling 
protects  components  from  high  gas  temperatures  by  blanketing 
the  component  with  a  film  of  cooler  air.  This  film  insu- 
lates the  component  and  acts  as  a  heat  sink.  The  problem 
encountered  with  film  cooling  is  the  prediction  of  the 
component's  wall  temperature.  Two  different  convective  heat 
transfer  coefficients  are  currently  used  for  this 
prediction: 

q"  =  h^(T    -  T  )  (1) 


f   aw    w 


and 


q"  =  hfT    -  T  )  (2) 

Here,  different  temperature  differences  are  used  to  define  h 
and  h£. 

The  objective  of  this  thesis  is  to  find  a  relation 
between  the  two  heat  transfer  coefficients  when  the  effects 
of  variable  properties  are  important. 
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II.  BACKGROUND 

A.   TURBINE  BLADE  COOLING  FUNDMENTALS 

For  a  given  size  gas  turbine,  raising  the  turbine  entry 
temperature  results  in  a  higher  specific  power.  This 
enables  the  designer  to  operate  the  'engine  cycle'  at  a 
higher  pressure  ratio,  resulting  in  improved  thermal  effi- 
ciency of  the  thermodynamic  cycle  and  improved  overall  effi- 
ciency. The  net  result  is  a  lower  specific  fuel 
consumption. 

Raising  the  turbine  entry  temperature  results  in  higher 
component  temperatures.  These  components  become  more  prone 
to  failure  at  higher  temperatures  through  greater  oxidation 
rates  of  their  surfaces,  excessive  creep,  and/or  by  thermal 
fatigue.  Figure  1  (Hennecke,  1982)  shows  a  comparsion  of 
the  rising  turbine  entry  temperature  to  the  rising  blade 
material  temperature  limitations.  Notice  that  uncooled 
metal  components  are  limited  to  a  turbine  entry  temperature 
of  approximately  1200  K.  Current  fuel  economy  goals  require 
a  turbine  entry  temperature  of  1600  K  and  greater. 
Numereous  cooling  systems  have  been  devised  to  maintain 
internal  components  at  or  below  their  maximum  allowable 
temperature.  Vanes,  blades,  disks,  liners,  and  casings  may 
all  require  some  form  of  cooling  protection. 

Figure  2  (Daniels  and  Schultz,  1982)  illustrates  two 
possible  cooling  medium  schemes.  Liquid  systems  use 
external  sprays  and/or  internal  convection  methods.  Air 
schemes  use  external  films  and/or  internal  convection 
methods.  Liquid  systems  have  several  thermodynamic  advan- 
tages over  air  systems.  Dispite  these  advantages,  however, 
air  systems  have  been  universally  accepted  due  to  their 
relative  mechanical  simplicity  to  implement.  Three  prin- 
ciple methods  '  of  air   cooling  currently   exist.    They   are 


11 


internal   convection,    film,    and   transpiration/effusion. 
Figure  3  (Hennecke,  1982)    illustrates  these  characteristic 
types  of  cooling  for  turbine  blades. 

Internal  convection  is  cooling  by  convection  between  the 
internal  sections  of  the  blade  and  the  air.  This  type  of 
cooling  may  be  subdivided  into  several  methods  consisting 
of: 

(1)  Straight  through  radial  holes  or  channels  of  a  single 
or  multipass  nature; 

(2)  Enhanced  heat  transfer  due  to  an  enlargement  of  the 
internal  surface  area,  promotion  of  . turbulance,  or 
impingement ; 

(3)  An  insert  that  distributes  the  air  according  to  the 
localized  cooling  requirements. 

Film   cooling   is   localized   ejection   of   air   through 

discrete  slots   or  holes  on  the   blade.    This  results   in  a 

film  of  air  around  the  blade  which  is  cooler  than  combustion 

gas  temperature,  acting  as  a  heat  sink.    The  film  also  acts 

like  a  thermal  barrier  between  the  hot  gases  and  the  blade's 

surface.     This  type   of   cooling  may  be   subdivided   into 

several  methods  consisting  of: 

(1)  Single  films;  where  specific  parts  of  the  blade  are 
covered  by  a  ^Im; 

(2)  Full  coverage;  where  many  films  are  placed  close 
enough  together  so  they  cover  part  or  all  of  the 
blade. 

Transpiration  or  effusion  cooling  is  the  bleeding  of  air 
through  the  blade  at  every  point   to  form  a  protective  film. 
The  blade   is  made  of  a   porous  material  to   accomplish  this 
method  of  cooling. 

Typically,  cooling  air  is  extracted  from  various  stages 
in  the  compressor  and  supplied  to  the  components  that 
require  cooling.  Figure  4  (Hennecke,  1982)  illustrates  the 
coolant  flow  path  for  a  single  stage  of  a  typical  gas 
turbine.  Typical  coolant  flow  rates  as  a  percentage  of  the 
gas  turbine's  total  flow  rate  are  included.  Notice  that  the 
guide  vanes  and  turbine  blades  have  the  highest  coolant  flow 
rates.  '  These  two  components  are  of  primary  concern  for  they 
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have  the  largest  amount  of  surface  area  exposed  to  the  main- 
stream flow  and  consequently  the  highest  thermal  loading. 

Current  applications  use  a  combination  of  the  internal 
convection  and  film  cooling  methods.  Figure  3  (Daniels  and 
Schultz ,  1982)  illustrates  a  typical  guide  vane  cooling 
configuration.  Cooling  air  comes  in  from  the  top  and  bottom 
of  the  vane.   Notice  the  internal  convection  methods  used: 

(1)  Air  is  distributed  to  the  internal  surface; 

(2)  Impingement  on  the  surface  after  distribution; 

(3)  Ribs  and  pedestals  used  to  enlarge  the  surface  area. 
After  the  air  has  cooled  the  vane  by  internal  convection,  it 
is  ejected  into  the  mainstream  and  used  as  a  film  on  the 
exterior  surface.  Initial  stages  utilize  the  internal 
convection  and  film  cooling  principles  as  illustrated  in 
Figure  5.  Latter  stages  may  employ  internal  convection 
cooling  solely  due  to  the  lower  gas  temperatures  encountered 
in  these  stages. 

The  goal  in  using  any  cooling  system  is  to  maintain  the 
internal  components  temperatures  at  or  below  their  maximum 
allowable  temperature.  Achievement  of  this  goal  will  result 
in  high  reliability  and  long  service  life.  In  addition,  it 
is  desirable  to  obtain  (as  closely  as  possible)  a  uniform 
blade  temperature  distribution  to  minimize  thermal  stresses. 

In  the  selection  of  a  cooling  system  to  be  employed,  two 
key  limitations  must  be  kept  in  mind.  These  limitations 
are: 

(1)  The  minimal  ^  amount  of  cooling  air  should  be  used  to 
minimize  it's  detrimential  efrect  on  engine  perform- 
ance ; 

(2)  The  effectiveness  of  a  particular  cooling  method  to 
keep  the  component  sufficiently  cooled'  for  the 
desired  turbine  entry  temperature. 

The  detrimental  effect  on  engine  performance  is  two 
fold.  First,  extracting  air  from  the  compressor  reduces  the 
available  output  power.  Some  work  can  be  recovered  by  rein- 
troducing  the  cooling   air  back   into  . the  mainstream   flow 
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after  cooling.  Second,  ejection  of  air  through  the  blade 
has  a  detrimential  effect  on  the  aerodynamic  efficiency  of 
the  turbine.  Diminished  aerodynamic  efficiency  results  from 
pressure  losses  and  flow  separation  m  the  wakes  trailing 
the  blades. 

Figure  6  (Hennecke,  1982)  shows  the  effectiveness  of 
different  methods  used  to  cool  components  as  a  function  of 
the  turbine  entry  temperature.  'Cooling  effectiveness'  is 
defined  as  the  ratio  of  actual  to  theoretical  blade  tempera- 
ture reduction.  This  'cooling  effectiveness'  is  not  only  a 
function  of  the  geometric  configuration  involved,  but  also 
the  coolant  flow  rate.  As  the  turbine  entry  temperature 
increases,  so  must  the  required  cooling  effectiveness. 

Figure  7  shows  the  cooling  flow  required  by  various 
cooling  methods  as  a  function  of  the  turbine  entry  tempera- 
ture. Notice  that  convection  cooling  requires  the  most 
cooling  flow  rate  while  effusion  requires  the  least  at  any 
given  turbine  entry  temperature.  Although  this  figure  is 
for  specific  leading  edge  and  coolant  temperatures,  the 
trend  is  the.  same  for  other  temperature  combinations. 

Figures  6  and  7  indicate  that  to  achieve  the  highest 
possible  turbine  entry  temperature  with  minimum  detrimential 
effects,  effusion  cooling  should  be  used.  Effusion  cooling, 
however,  is  not  currently  used  due  to  structural  and  manu- 
facturing difficulties,  and  because  the  porous  blades  often 
become  clogged.  This  is  one  reason  why  film  cooling  methods 
are  currently  employed. 

Logically,  one  could  assume  that  the  more  coolant  pumped 
into  the  film  and  the  fuller  the  coverage,  the  higher  the 
turbine  entry  temperature.   However,  there  are  limitations: 

(1)  The  detrimential  aerodynamic  effects  of  film  cooling 
may  offset  the  gams  of  a  higher  turbine  entry 
temperature ; 

(2)  "Lift-Off"  may  occur  where  the  film  ceases  to  protect 
the  surface. 
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Thus,   thermal  design  of  vanes   and  blades  with  film  cooling 
must  be  optimized. 

B.   FILM  COOLING 

Film  cooling  is  the  ejection  of  coolant  into  the 
boundary  layer  on  the  blade.  Figure  8  shows  a  typical  film 
cooling  geometry  with  temperature  profiles  for  a  row  of 
holes.  Notice  on  the  top  view,  at  some  distance  downstream, 
the  flows  from  adjacent  jets  converge  to  form  a  continuous 
film  over  the  blade's  surface.  In  this  region,  lateral 
variations  become  negligible,  allowing  a  two  dimensional 
analysis  to  be  conducted  with  reasonable  accuracy. 

The  typical  temperature  profiles  on  Figure  8  illustrate 
the  effect  of  film  cooling  on  heat  transfer.  The  left 
profile  is  without  film  cooling.  Heat  transfer  is  from  the 
freestream  to  the  blade's  surface.  The  driving  force  is  a 
function  of  the  freestream  recovery  temperature  and  wall 
temperature.  The  right  profile  is  with  film  cooling.  Heat 
transfer  is  from  the  freestream  to  the  coolant  and  from  the 
coolant  to  the  blade's  surface.  The  driving  force  is  a 
function  of  the  freestream  recovery  temperature,  coolant 
recovery  temperature,  and  wall  temperature.  As  the.  flow 
moves  downstream  from  the  injection  site,  the  coolant  mixes 
with  the  existing  boundary  layer.  Eventually  a  temperature 
profile  similar  to  the  left  will  prevail.  When  this 
happens,  the  heat  sink  and  insulating  effects  of  the  film 
become  minimal. 

Figure  9  illustrates  a  typical  heat  transfer  coefficient 
and  adiabatic  wall  temperature  distribution  around  a  turbine 
blade.  Notice  that  the  highest  temperature  and  heat 
transfer  coefficient  occurs  at  the  nose  of  the  blade,  near 
the  stagnation  point,  and  the  location  of  the  highest 
thermal  loading.  Because  of  this,  injection  sites  are 
located  in  the  nose  and  utilize  the  lowest  coolant  tempera- 
ture.   In  the  region  around  the  stagnation  point,   the  heat 
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transfer  coefficients  in  Figure  9  are  representative  of  a 
laminar  boundary  layer.  At  approximately  one- fourth  the 
distance  along  the  suction  surface,  transition  to  a  turbu- 
lent boundary  layer  begins,  resulting  in  an  increase  in  the 
adiabatic  wall  temperature  and  heat  transfer  coefficient 
approximately  three-fourths  the  distance  along  the  suction 
surface.  The  heat  transfer  coefficients  along  the  pressure 
side  in  Figure  9  are  representative  of  a  turbulent  boundary 
layer  developing  with  a  constant  streamwise  pressure 
gradient  (after  an  initial  drop  from  the  stagnation  point). 
Near  the  mid-chord  position,  the  heat  transfer  coefficients 
begin  to  rise  as  a  strong  favorable  pressure  gradient 
develops. 

C.   BOUNDARY  LAYER  EQUATIONS 

To  predict  the  heat  flux  from  the  mainstream  to  the 
blade,  the  equations  governing  the  physical  situation  must 
be  solved.  Figure  10  illustrates  the  coordinate  system 
used.  Taking  a  differential  element  and  applying  the 
conservation  principles  of  mass,  momentum,  and  energy;  the 
governing  equations  are  derived.  Using  the  governing  equa- 
tions, a  solution  to  the  heat  transfer  problem  is  consider- 
ably difficult  since  flow  and  temperature  distributions  are 
very  complex.  These  equations  are  signif iciently  simplified 
by  making  the  following  assumptions: 
(  1)   u  >>  V 

(2)  —  -  0 

3T      3T 

(3)  >>>  

dy  8x 


(4) 


2 

X  +  -u  =  0 
3 


(5)  No  internal  heat  sources  exist 

(6)  No  chemical  reactions  occur  within  the  element 
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(7)  Body  forces  are  negigible 

(8)  No  work  is  performed  by  external  fields 

Applying  these  assumptions  reduces  the  governing  equa- 
tions to  the  boundary  layer  equations.  For  turbulent 
boundary  layers,  the  following  form  is  for  time-averaged 
velocities  and  variable  properties: 


CONTINUITY 


MOMENTUM 


3  ___    a  _ 

(pu)  * (pv)  =0  C3) 

3x        3y 


3u        3u     dP    3    3u 


(pu)—  +  (pv)—  = *  —iu pu'v')      (4) 

3x       3v     dx    3v   3y 


ENERGY 


_^  3H    3H   _dP    3q     STT 

(pu) —  +  (pv) —  =  u —  +  —  +  T —  (5) 

3x    '   3y     dx    3y     3y 


where 


dU 


and 


T  =  u— -  PU  '  V  '  (6) 

9y 


3T 


q  =  k pv'H'  (7) 

3y 


The  boundary  layer  equations,   as  presented,   are  highly 
nonlinear.    The   momentum  and  energy  equations   are  coupled 
through  the  temperature  dependency  of  density  and  viscosity. 
Thermal  conductivity   is  also   temperature  dependent.     The 


17 


dependence  of  these  transport  properties  on  temperature 
results  in  different  velocity/ temperature  profiles  and 
different  heat  transfer  coefficients  from  those  obtained  by- 
assuming  properties  which  are  constant  through  the  boundary 
layer. 

Assuming  constant  properties  exist   through  the  boundary 
layer,  the  boundary  layer  equations  may  be  rewritten  as: 


CONTINUITY 


9u    3v 

—  *  ~  =  0  (S) 

3x    3y 


MOMENTUM 


ENERGY 


au     3u      IdP    3    3u    

^_  ^    ^. —  = _  +  (^ u'v')  (9) 

3x     3y      pdx  •   3y   3y 


3T     3T    3    3T 


u— ^+  V —  =  —(a—  -  v'T')  CIO) 

3x     3y    3y   3y     • 


The  momentum  and  energy  equations  ((9)  and  (10))  are  uncou- 
pled. This  means  that  the  momentum  equation  can  be  solved 
independently  of  the  energy  equation.  In  addition,  the 
energy  equation  is  linear.  These  two  items  simplify  the 
process  of  solving  the  boundary  layer  equations. 

A  number  of  existing  methods  to  account  for  variable 
properties  in  boundary  layers  without  film  cooling  will  now 
be  discussed. 

Mainstream  gases  consist  of  air  and  the  combustion 
process  by-products  (principally  Carbon  Dioxide,  Carbon 
Monoxide,  and  water  vapor).  Viscosity  and  thermal  conduc- 
tivity of  this  mixture  can  be  determined  though  the  semiem- 
pirical  formula  from  White  (1974)  : 
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n   m-u- 

y  .   =  z LJ: (in 

mix   .  ,  ,    ^ 

1  =  1  n 

I   m  .  4) .  . 


j  =  l 


where 


(1  *   (u^/u^)^''^(w./w,)^/''} 


i 
This  formula  is  for  a  mixture  of  n  gases.  W-  are  the  molec- 
ular weights  and  m^  are  the  mole  fractions.  The  exact  same 
formula  is  recommended  for  thermal  conductivity  (y  is 
replaced  by  k).  Figure  11  (White  (1974))  illustrates  the 
temperature  variations  of  viscosity  and  thermal  conductivity 
for  these  gases  at  low  pressures.  Notice  that  for  a  200  R 
temperature  variation,  the  variations  in  viscosity  and 
thermal  conductivity  are  not  very  dramatic.  A  power  law  of 
the  form  (Kays  and  Crawford  (1980)): 

y  .      k  .  -     \  T  / 

mix     mix       e 
e       e 

can  be  used  to  determine  their  temperature  variations.    For 

small  temperature  differences  through  boundary  layers,   say 

20  C  or  less,  viscosity  and  thermal  conductivity  may  each  be 

considered  to  be  constant. 

By  considering  the  mainstream  gases   to  be  a  perfect  gas 

and  from  uniform  static  pressure  through  the  boundary  layer, 

the  density  variation  becomes   inversely  proportional  to  the 

absolute  temperature. 

P     T, 


-  (14) 


^e    T 
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Comparsion  of  equations  13  and  14  shows  that  temperature 
variations  result  in  larger  density  ratio  changes  than  \i  or 
k  changes. 
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III.  PRELIMINARY  CONSIDERATIONS 

The  objectives  in  the  thermal  design  of  gas'  turbine 
blades  are  to  obtain  an  acceptable  stress  level  and  design 
life.  To  do  this,  a  detailed  blade  temperature  distribution 
must  be  obtained. 

The  thermal  design  of  a  cooled  turbine  blade  is  a 
complicated  energy  balance.  The  design  process  can  be 
broken  down  into  three  major  categories  outlined  below  (in 
sequential  order): 

(1)  Predicting   the  heat   flux   from   the  hot   mainstream 
gases  to  the  blade's  surface; 

(2)  Performing   a  conducljiion   analysis   for  the   detailed 
mapping  or  the  blade  s  temperature  distribution; 

(3)  Predicting   the   heat   flux  from   the   blade   to   the 
internal  cooling  air. 

This  thesis  only  addresses  the   first  category:    a  relation 

is   developed  and   verified  between   the  two   types  of   heat 

transfer  coefficient  used  to  predict  wall  temperature. 

A.   CONSTANT  PROPERTY  CONSIDERATIONS 

Equations  1  and  2  are  used  to  calculate  the  blade's  wall 
temperature  after  the  heat  flux  is  obtained  by  solving  the 
boundary  layer  equations. 

q"  =  h^(T    -  T  )  CD 

^      f   aw    w^  ^    ■' 

and 

q"  =  hCT    -  T  )  (Z) 

*        o«     w 

Both  equations  reflect  the  presence  of  film  cooling. 

In  developing  thermal  boundary  layers  with  film  cooling, 
an  adiabatic  wall  boundary  condition  can  be  used  to  describe 
the  convection  process.  Using  this  concept  reduces  the 
driving  potential  from  a  three  temperature  potential  problem 
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to  a  two  temperature  potential  problem.  Equation  1  repre- 
sents this  concept.  Here,  hr  is  an  isoenergtic  heat 
transfer  coefficient.  The  adiabatic  wall  temperature,  T  , 
is  the  wall  temperature  when  there  is  no  -heat  transfer  to 
the  wall.  Adiabatic  wall  temperature  is  usually  expressed 
nondimensianally  in  terms  of  the  adiabatic  film  cooling 
effectiveness,  n^ 

'aw     r°> 

n  =  CIS) 

'rc    'r~ 

Since  the  entire  surface  is   considered  adiabatic,   T    is  a 

aw 

global  parameter.  As  a  result,  equation  1  is  dependent  on 
the  entire  convection  process  and  the  thermal  boundary 
layer. 

Compared  to  the  large  temperature  differences  between 
the  mainstream  gases  and  the  blade's  surface,  the  surface 
temperature  distribution  can  be  considered  relatively  isoth- 
ermal. An  isothermal  wall  is  ideal  since  thermal  stresses 
are  minimal.  Equation  2  represents  this  concept.  Here,  h 
is  defined  by  the  freestream  stagnation  temperature  and  the 
wall  temperature  difference  for  a  specific  coolant  tempera- 
ture. The  coolant  temperature  is  usually  expressed  nondi- 
mensionally  as  9: 


(1(3) 


There  has  been  considerable  controversy  over  these  two 
methods  of  predicting  the  wall's  temperature.  Two  key 
issues  of  this  controversy  deal  with  how  the  temperature 
difference  is  defined  and  the  effects  of  changing  coolant 
parameters  (temperature,  density,  mass  flow  rate,  etc.).  In 
equation  1,  when  the  temperature  difference  equals  zero, 
there  is  no  heat  transfer.    This  is  not  necessarily  true  if 
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T 
re 

-     T 

T      - 
w 

T 

equation  2  is  used.  To  account  for  viscous  dissipation,  the 
freestream  stagnation  temperature,  Tq^o,  is  replaced  by  the 
freestream  recovery  temperature,  "^-j-oo-  ^^"*®  heat  transfer 
coefficient,  h,  is  then  redefined  as: 

q"  =  h(T    -  T  )  [1"^ 


Changing  the  film  cooling  parameters  will  have  different 
effects  on  these  equations.  Generally,  the  recovery  temper- 
ature is  determined  from  freestream  conditions  (temperature, 
velocity,  density),  however,  if  it  is  based  on  near  wall 
conditions,  the  recovery  temperature  in  a  film  cooled 
boundary  layer  is  the  same  as  the  adiabatic  wall  tempera- 
ture. In  equation  1,  both  h^  and  T  are  affected,  whereas, 
in  equation  17  only  h  is  affected. 

For  a  given  location,  the  heat  flux  will  be  the  same  for 
either  method.  As  such,  equation  1  is  set  equal  to  equation 
17  and  a  relation  between  the  heat  transfer  coefficients 
exists.  Using  the  definitions  of  ri  and  9,  and  dividing  by 
the  heat  transfer  coefficient  for  the  case  of  no  film 
cooling,  then  yields: 

h     h^ 

—  =  —(1  -  n9)  (18) 

o     o 

This  equation  relates  adiabatic  and  isothermal  wall  boundary 
conditions  for  properties  that  are  constant  through  the  film 
cooled  boundary  layer.  This  relationship  is  plotted  on 
Figure  12  (Ville,  Cunat ,  and  Richards,  1978).  Note  the 
linearity  between  h/h^  and  9.  This  linearity  is  due  to  the 
fact  that  the  energy  equation  is  linear  for  constant  proper- 
ties. This  equation  is  valid  for  a  given  dov;nstream  loca- 
tion, where  adiabatic  conditions  are  represented  by  the 
horizontal  intercept,  and  isothermal  conditions  by  the  line. 
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On  Figure  12-,  the  axis  intercept  points  have  special 
meaning.  The  h/hQ  axis  intercept  point  represents  the  heat 
transfer  coefficient  when  the  coolant  is  injected  at  the 
same  temperature  as  the  mainstream  (T^|-=Tpco»  6=^0).  The  9 
axis  intercept  point  indicates  a  wall  condition  that  is 
locally  adiabatic.  At  this  point,  9  equals  the  inverse  of 
the  effectiveness,  n. 

Figure  12  is  only  good  for  one  downstream  location.  As 
you  move  further  downstream,  the  slope  of  the  curve  becomes 
less  and  may  eventually  become  horziontal  indicating  no 
cooling  effect  from  the  film.  This  is  a  result  of  a  rising 
film  temperature  due  to  heat  transfer,  diffusion,  and  mixing 
with  mainstream  gases.  As  the  'film  temperature  rises,  the 
value  of  9  decreases  and  the  value  of  n  increases. 

B.   VARIABLE  PROPERTY  CONSIDEFJ^TIOXS 

With  variable  properties,  the  momentum  and  energy 
boundary  layer  equations  become  coupled  because  of  the 
density  term.  In  addition,  the  energy  equation  becomes 
highly  nonlinear.  Because  of  this,  the  linearity  exhibited 
on  Figure  12  no  longer  exists.  Ligrani  and  Camci  (1985) 
illustrated  this  nonlinearity  that  exists  between  h/h^  and 
9.  They  also  derived  an  empirical  relation  to  relate  the 
constant  property  data  of  Jabbari  and  Goldstein  (1978)  to 
their  variable  property  data. 

To  find  the  relation  between  h/h  and  9  for  the  variable 
property  case,  we  start  with  the  defining  equations  for  the 
heat  transfer  coefficients. 

For  variable  properties  through  the  boundary  layer,  the 
variable  property  isoenergetic  heat  transfer  coeffi- 
cient, h ^,  is  defined  using: 

q"  =  hj(T^,  -  TJ  (19) 

The  heat  transfer  coefficient,  h,  is  then  defined  using: 

q"  =  h(T^^  -  T^)  (20) 
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The  heat  flux  is  now  obtained  by  solving  the  variable  prop- 
erty boundary  layer   equations.    At  a  given   location,   the 
heat   flux   from  equations   19   and   20   will  be   the   same. 
Equating  these  gives: 

hr(T^,  -  T,,)  =  h(.T    -  T  )  (21) 

f   aw     w^     ^  r°o     w 

Now  the  adiabatic  wall  temperature   is  expressed  in  terms  of 
an  adiabatic  film  cooling  effectiveness  defined  as: 

1  .   T 

n  =  (22) 

T     -  T 

'rc    "r°° 

9  is  defined  as  before   (equation  16).    By  substitution  and 
rearrangement  the  variable  property   relation  between  h/h 
and  9  then  becomes: 

h      h 

—  =  -i(l  -  ne)  (25) 

O  0 

This  equation  is  for  one  downstream  location  at  one 
value  of  the  density  ratio,  P^^/Poo*  ^^*^  °^^  value  of  any  of 
the  following ;. blowing  ratio, M;  momentum  flux  ratio,  I;  or 
the  velocity  ratio,  u_/u^.  A  different  line  results  from 
equation  23  for  each  value  of  the  effectiveness,  n ,  and  the 
normalized  heat  transfer  coefficient,  h£/hQ.  The  effective- 
ness, n,  may  have  a  strong  dependence  on  the  density  ratio, 
Pp/ Poo»  ^^  indicated  by  the  single  row  film  cooling  data  of 
Pedersen,  Eckect,  and  Goldstein  (1977). 

Figure  13  is  a  domain  map  showing  density  ratio  verses 
blowing  ratio  of  the  Pedersen,  Eckect,  and  Goldstein  (1977) 
data.  The  box  symbols  indicate  that  effectiveness  has 
little  or  no  dependence  on  the  density  ratio  for  low  blowing 
ratios.  The  circle  symbols  indicate  that  effectiveness  has 
a  strong  dependence  on  the  density  ratio.  The  triangles 
indicate  that  the  effect  of  the  density  ratio  on  the  effec- 
tiveness and  heat  transfer  is  highly  nonlinear. 
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Figure  14  shows  one  possible  effect  of  variable  proper- 
ties on  h/h  verses  8  data.  Here,  each  of  the  three  linear 
curves  represents  equation  23  with  a  different  t\  .  The  same 
h^/h  was. used  for  each  line  to  simplify  Figure  14.  Now, 
each  line  is  also  for  one  value  of  the  density  ratio,  P(- /poo , 
which  is  equivalent  to  T^/Tq.  The  variation  of  the  density 
ratio  with  6  is  then  given  by: 

=  9  (241 

T  /T    -  1 
w   r°o 

if   T  /T„  =  T   /T    for  constant   T  '.         The   9   for  a   given 

c '  *     re   r<»  w  '^ 

density  ratio  will  then  be  a  straight  vertical  line 
eminating  from  the  horizontal  axis.  The  intersection  of 
this  vertical  line  and  equation  23  will  be  where  a  given 
h./h  data  point  is  located.  The  curved  line  represents  the 
collection  of  the  h/h  verses  0  data  points.  It  curves 
downward  because  of  the  way  the  data  collects.  As  the 
density  ratio  increases  (vertical  lines  move  to  the  right), 
the  effectiveness  increases.  This  causes  the  9  axis  inter- 
cept point  (l/n)  to  move  to  the  left. 
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IV.  RESULTS 

Using  equation  23,  the  adiabatic  effectiveness  was  used 
to  deduce  1^/h  using  the  data  of  Pedersen,  Eckect,  and 
Goldstein  (1977)  for  one  row  of  injection  holes.  In  addi- 
tion,  the  adiabatic  effectiveness  was  deduced  from  h/h  and 
hr/h  from  the  data  of  Ligrani  and  Camci  (1985)  for  two  rows 
of  injection  holes. 

In  the  following  sections,  the  results  are  presented. 

A.   DATA  OF  PEDERSEN,  ECKECT,  AND  GOLDSTEIN  (1977) 

Figure  15   shows  the   Pedersen,   Eckect,    and  Goldstein 

(1977)    data  used   to  deduce   h/h^.    Figure   13  gives   the 

blowing  ratio  and   the  density  ratio  values  for  which  n  was 

selected  for  deducing  h/h  . 

A  value  of  1.  0  was  assumed   for  the  h/h   axis  intercept. 

'  o  ^ 

Although  in  many  cases,  a  reasonable  assumption,  Metzger, 
Carper,  and  Swank  (1968)  showed  that  this  axis  intercept  may 
not  be  equal  to  1.  0.  They  showed  that  the  injection  has  an 
effect  on  the  heat  transfer  coefficient  even  when  the 
coolant  is  injected  at  the  freestream  temperature,  primarily 
due  to  boundary  layer  thicking. 

Normalized  heat  transfer  coefficients  as  a  function  of  9 
are  shown  on  Figures  16  thru  19  for  the  Pedersen,  Eckect, 
and  Goldstein  (1977)  data.  Four  axes  are  used  to  denote 
four  downstream  locations  (x/d).  The  light  linear  lines 
represent  equation  23  for  a  different  value  of  effective- 
ness, T].  The  dark  lines  represent  the  collection  of  h/h^ 
verses  9  data  points. 

Figure  16  is  for  a  blowing  ratio,  M,  of  0.213.  This 
figure  shows  that  h/h^  verses  9  is  steeper  at  small  x/d  and 
becomes  flatter  as  the  distance  from  the  injection  site 
increases  (increasing  x/d).  This  figure  shows  reasonable 
linearity  with  some  nonlinearity  near  the  injection  site. 
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Figure  17  is  for  a  blowing  ratio,  M,  of  0.515.  This 
figure  shows  the  same  trends  as  with  Figure  16.  Only  now, 
the  nonlinearity  near  the  injection  site  is  more  pronounced 
and  is  also  evident  at  a  location  further  downstream. 

Figure  18  is  for  a  blowing  ratio,  M,  of  1.05.  This 
figure  shows  that  the  variation  of  h/h^^  with  9  is  nonlinear. 
As  the  distance  from  the  injection  site  increases 
(increasing  x/d) ,  the  h/h  variation  with  9  becomes  more 
horizontal.  Notice  that,  in  general,  at  a  particular 
blowing  ratio,  the  slope  becomes  more  negative  as  9 
increases. 

Figure  19  is  for  a  blowing  ratio,  M,  of  1.96.  This 
figure  shows  a  highly  nonlinear  behavior  due  to  the  'lift- 
off of  film  cooling  jets  from  the  wall,  and  possible  pene- 
tration through  the  boundary  layer. 

B.   DATA  OF  LIGRANI  AND  CAMCI  (1985) 

Figure  20  (Ligrani  and  Camci,  1985)  illustrates  the 
downward  sloping  trend  for  variable  property  h/h^  verses  9 
data.  This  figure  also  shows  the  linear  h/h-  verses  9  data 
(constant  property  case).  Notice  the  similarity  between  the 
variable  property  trend  here  and  the  trends  of  the  Pedersen, 
Eckect,  and  Goldstein  (1977)  data  on  Figure  18. 

Using  h/h  and  h^/h  from  Ligrani  and  Camci  (1985),  n 
was  deduced  from  equation  23.  Results  are  plotted  as  a 
function  of  x/d  in  Figures  21  thru  31,  and  as  a  function 
of  ^  in  Figures  32  thru  42.  The  later  function,  C,  is  a 
downstream  distance  normalized  by  the  blowing  ratio.  It  is 
defined  by: 

X  +  1  .909.1      u^    ^  .. 

MS        ^  u 

JO 

where  S=7id/6  is   the  width  of  an   equivalent  two  dimensional 
slot  (Jabbari  and  Goldstein,  1978)  .   Goldstein  (1971)  gives 
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two  equations  to  -predict  laterally  averaged  effective- 
ness,  n>  for  two  dimensional  film  cooling.  These 
predictions  are  given  by; 

n  =  (1  +  0.2490"°-^  (26) 


and 


l.DPr^^^ 


n  =  no 

1  +  0.329C   e 


(27) 


where 

6  =  1*  I-SCIO""^)  (Re^(-S.))sinY  (28) 

Woo 

These  predictions   are   plotted  on  the  n  verses  ^  graphs 
(Figures  32  thru  42). 

Figures  21  thru  42  are  for  different  blowing  ratios  (M  = 
1.242,  1.003,  0.481,  0.661,  1.557),  momentum  flux  ratios  (I 
=  1.382,  0.925,  0.417),  and  velocity  ratios  (u  /u^  =  0.630, 
1.113,  0.922).  In  general,  these  figures  show  that  as  the 
blowing  ratio,  momentum  flux  ratio,  and  velocity  ratio 
increase,  the  adiabatic  film  cooling  effectiveness 
increases.  The  largest  increase  in  n  occurs  for  the  highest 
value  of  0.  This  corresponds  to  the  case  of  the  highest 
density  ration.  On.  Figures  21  thru  31,  the  crossover  point 
between  the  01  and  02  curves  moves  downstream  (from  x/d  =  20 
to  x/d  =40)  as  the  momentum  flux  ratio  and  velocity  ratio 
increase.  On  Figures  32  thru  42,  as  the  blowing  ratio, 
momentum  flux  ratio,  and  the  velocity  ratio  decrease,  the 
data  moves  to  a  higher  value  and  spreads  out.  All  data 
points  on  the  n  verses  C  graphs  are  below  the  predictions  of 
equations  26  and  27  but  do  follow  the  same  trend  as  the 
predictions. 
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V.  SUMMARY 

A  relation  was  derived  that  relates  the  normalized  heat 
transfer  coefficients  h  and  h^.  These  heat  transfer  coeffi- 
cients are  for  isothermal  wall  boundary  conditions,  h,  and 
adiabatic  wall  boundary  conditions.  h£ .  Each  are  for  when 
variable  properties  are  considered  to  exist  through  the 
boundary  layer.   This  relation  was  presented  as: 

h       hr         /s 

—  =  -^(1  -  ne)  (25) 

o    o 

When  the   effects  of  variable  properties   are  important   in 

film  cooled  turbulent  boundary  layers,  the  variable  property 

effectiveness,   n ,   and  the  above  equation  are  valid  for  one 

density  ratio  and  blowing  ratio  only. 

The  normalized  heat  transfer   coefficient,   h/hQ ,    was 

deduced  from  the  Pedersen,  Eckect,  and  Goldstein  (1977)  data 

on  injection  with  one   row  of  holes.    The   adiabatic  film 

cooling  effectiveness,   n,   was  deduced  from  the  Ligrani  and 

Camci  (1985)  data  on  injection  with  two  rows  of  holes.    The 

data  deduced  from  equation  23   showed  the  nonlinear  trend  of 

h/h^  verses  6  and    as  the  density  ratio,  P_/p_,  varied, 
o  ^  ^ 

The  effect  of  density  ratio  on  h/ho  verses  9  results  in 
a  change  of  the  characteristic  curve  from  linear  to  nonli- 
near. This  change  is  not  always  evident.  It's  existance  or 
extent  is  dependent  upon  the  injection  parameters  (blowing 
ratio,  momentum  flux  ratio,  and  velocity  ratio)  and  the 
downstream  distance.  This  effect  is  due  to  the  nonlinearity 
introduced  into  the  energy  equation  and  the  coupliing  of  the 
momentum  and  energy  equations  when  variable  properties 
through  the  boundary  layer  are  considered. 
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VI.  RECOMMENDATIONS 

Experiments  need  to  be  conducted  to  determine  the  accu- 
racy and  dependability  of  the  proposed  variable  property 
relation  on  injection  system  geometry,  injection  flow  rate 
parameters,  fluid  mechanics  of  the  flow,  the  thermal 
boundary  layer  conditions,  and  other  parameters.  In  order 
to  completely  verify  the  proposed  equation,  simultaneous 
measurements  of  variable  property  effectiveness  and  variable 
property  heat  transfer  coefficients  (based  on  the  wall  and 
recovery  freestream  temperature  difference)  need  to  be 
conducted. 
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1.  Internal  Cooling 

-  straight  through 
or  multipass 

-  enhanced  heat  transfer 

-  insert 


2.  Film  Cooling 

-  single  films 

-  full  coverage 


<yi^ 


3.  Effusion  Cooling 


Fig.  3    Characteristic  Types  of  Blade 
Cooling  Configurations. 
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Fig.  4    Coolant  Flow  in  a  Typical 
High  Pressure  Turbine  Stage. 
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Fig.     5         Coolant    Flow   in   a   Typical 
Nozzle    Guide   Vane. 
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Fig.  6    Turbine  Entry  Temperature  Verses 
Required  Cooling  Effectiveness. 
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Fig.  7    Comparative  Cooling  Requirements, 

Leading  Edge  Temperature  1000  C. 

Cooling  Air  Temperature  260  C. 
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Main  Flow 


-^-^^   Injection  of  Cooling  Air 


Region  of  Single  Jets       Mixed  Region 
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Fig.  8    Film  Cooling  by  a  Row  of  Holes 
with  Sample  Temperature  Profiles. 
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Fig.     10         Curvilinear    System   of    Coordinates, 
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Fig.  11    Viscosity  and  Thermal  Conductivity 
of  Common  Gases  at  Low  Pressures. 
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Fig.  12    Linear  Variation  of  Heat  Transfer 

Coefficient  with  Nondimensional 

Coolant  Temperature. 
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Fig.  13    Domain  Map  of  Blowing  Ratio  and  Density 
Ratio  for  Pedersen,  Eckect,  and  Goldstein  (1977)  Data. 
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—  from  equation   4 
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from  equation   15 


Fig.  20    Variation  of  h/h  with  0  for  Two  Rows 
of  Injection  from  Ligrani  and  Camci  (1985). 
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